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| Abstract 

| We study the generalized Chaplygin gas model (GCGM) using Gamma-ray bursts as 

cosmological probes. In order to avoid the so-called circularity problem we use cosmology- 
independent data set and Bayesian statistics to impose constraints on the model parameters. 

■ We observe that a negative value for the parameter a is favoured in a flat Universe and the 
| estimated value of the parameter Hq is lower than that found in literature. 

PACS number: 98.80.Es, 98.70.Rz 

\ 1 Introduction 

One of the most important problems of Modern Cosmology is the determination of the matter 
content of the Universe. The rotation curve of spiral galaxies [I], the dynamics of galaxy clusters 
[2] and structure formation [3], indicate that there is about ten times more pressureless matter 
in the Universe than can be afforded by the baryonic matter. The nature of this dark matter 
component remains unknown. Moreover, the Type la supernovae (SNe la) data indicates that 
the Universe is accelerating [lj . Models considering matter content dominated by an exotic fluid 
whose pressure is negative [5], modified gravity theories such as f(R) [6] and the evolution of 
\ an inhomogeneous Universe model described in terms of spatially averaged scalar variables with 

matter and backreaction source terms [7] are some of the proposals to explain this current phase 

■ of the Universe. At the same time, the position of the first acoustic peak in the spectrum of CMB 
\ anisotropies, as obtained by WMAP, favours a spatially flat Universe [8]. Combining all these 

data and if we consider the matter content of the Universe dominated by a fluid with negative 
pressure we have a scenario with a proportion of Q m ~ 0.27 and fide ~ 0.73, with respect to the 
critical density, for the fractions of the pressureless matter and dark energy, respectively. This 
,£_J | scenario is usually called as the concordance cosmological model. 

^ The question is to know what is the nature of the dark matter and dark energy components. 

For dark matter many candidates have been suggested such as axions, a particle until now 
undetected which would be a relic of a phase where the grand unified theory was valid [9], 
the lightest supersymmetric particle (LSP) like neutralinos [10] and the Kaluza-Klein particles 
[llj that are stable viable Weakly Interacting Massive Particles (WIMPs) and arise in two 
frameworks: In Universal Extra Dimensions [12] and in some warped geometries like Randall- 
Sundrum [13J. For the dark energy, in the hydro dynamical representations of matter, the most 
natural candidate is a cosmological constant, but there is a discrepancy of some 120 orders of 
magnitude between its theoretical and observed values [14]. For this reason, other candidates 
have been suggested like quintessence models that involve canonical kinetic terms of the self- 
interacting scalar field with the sound speed <? s = 1 [15] and k-essence models that employ 
rather exotic scalar fields with non-canonical (non-linear) kinetic terms which typically lead to 
a negative pressure [16]. More recently, a string-inspired fluid has been evoked: The Chaplygin 
gas p__], that appears as a promising candidate for the dark sector of the Universe. 
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The Chaplygin gas is represented by the equation of state 

Pc = ~- , (1) 
Pc 

where p c represents the pressure, p c the fluid density and A is a parameter connected with the 
sound speed. This equation of state is suggested by a brane configuration in the context of 
string theories [18]. However, a more general equation of state has been suggested [l9j: 

Pc = ~— , 2 

Pc 

where again p c and p c stand for the generalized Chaplygin gas component and a is a new 
parameter, which takes the value 1 for the traditional Chaplygin gas but values larger than 1, 
or even negative may be considered. This is the so-called generalized Chaplygin gas. 

Much observational data that has been used for comparison with the theoretical cosmological 
models like the generalized Chaplygin gas model (GCGM). The spectra of anisotropy of cosmic 
microwave background radiation [20j, baryonic acoustic oscillations |21| . the integrated Sachs- 
Wolfe effect [22], the matter power spectrum [23], gravitational lenses [21], X-ray data [25] 
and ages estimates of high- z objects [26] have been used in this sense. Also, constraints from 
combined data sources have been obtained in [27J. Another tool used to make this comparison 
is the Hubble diagram, the plot of redshift z versus luminosity distance cLl = ^/C/AttJ 7 , where 
C is the luminosity (the energy per time produced by the source in its rest frame) and J- is the 
measured flux, i.e., the energy per time per area measured by a detector. Normally, the SNe 
la data are considered good standard candles and are used to construct the Hubble diagram, 
because their luminosity are well known [H[28]. In particular, constraints on the Generalized 
Chaplygin gas have been studied in [29] . These assumptions rest on a foundation of photometric 
and spectroscopic similarities between high- and low-redshift SNe la. But this discussion is not 
yet finished [30] . The other problem comes from the fact that there still does not exist SNe la 
data with z > 1.8. To know the properties and behavior of dark energy for high values of z we 
will have to wait for new data of the SNe la or to find other distance indicators. In this sense, to 
extend the comparison between observational data and theoretical models at very high redshift 
we propose to use Gamma-ray bursts (GRBs) due to the fact that they occur in the range of 
high z beyond the SNe data found today [31]. 

The GRBs are jets that release ~ 10 51 — 10 53 ergs or more for a few seconds and becomes, 
in this brief period of time, the most bright object in the Universe. They were discovered in 
the sixties by the Vela satellites in the "Outer Space Treaty" that monitored nuclear explosions 
in space [32]. Launched in 1991 The Burst and Transient Source Experiment on the Compton 
Gamma-Ray Observatory (BATSE on the Compton GRO) [33] observations concluded that the 
angular distribution of the GRBs on the sky is isotropic within statistical limits. This study 
ruled out the idea that the GRBs are galactic objects, but it is consistent with the bursts being 
extra-galactic sources at cosmological distances. More recently, the SWIFT mission (launched 
in 2004) has provided the most accurate GRB data, available in the Swift BAT Catalog. 

The search for a self-consistent method to use the GRBs in cosmological problems is intense 
and promising. But the possibility of using GRBs as standard candles is not a simple question. 
GRBs are known to have several light curves and spectral properties from which the luminosity 
of the burst can be calculated once calibrated, and these can make GRBs into standard candles. 
Just as with SNe la, the idea is to measure the luminosity indicators, deduce the source lumi- 
nosity, measure the observed flux and then use the inverse-square law to derive the luminosity 
distance. The difficulty arises when these indicators are a priori established through some cos- 
mological model like the concordance one. This means that the parameters of the calibrated 
relations of luminosity /energy are still coupled to the cosmological parameters derived from a 
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given cosmological model. This is the so called circularity problem. This problem appears in 
several works that have made use of these GRBs luminosity indicators as standard candles at 
very high redshift [M]. It is possible to treat the circularity problem with a statistical approach 
|35j . On the other hand, many papers have dealt with the use of so called Amati relation, or 
the Ghirlanda relation for this purpose [36]. However, as argued recently in [37], these pro- 
cedure involve many unjustified assumptions which if not true could invalidate the results. In 
particular, many evolutionary effects can affect the final outcome. However, recently Liang et 
al. [38\ [39j HQ] made a study considering SNe la as first-order standard candles for calibrating 
GRBs, the second-order standard candles. The sample in reference [38] was calibrated from the 
192 supernovae obtained in [41j . The updated sample used in [39], HQ] has been obtained and 
calibrated cosmology-independently from the Union2 (557 data points) compilation [32] released 
by the Supernova Cosmology Project Collaboration. In these articles the authors found rele- 
vant constraints on the Cardassian and Chaplygin gas model by adding to the GRB data the 
SNe la (Union2), the Shift parameter of the Cosmic Microwave Background radiation from the 
seven-year Wilkinson Microwave Anisotropy Probe and the baryonic acoustic oscillation from 
the spectroscopic Sloan Digital Sky Survey Data Release galaxy sample. The sample obtained 
in [39] will be used in our analysis. These authors obtain the distance moduli p of GRB in 
the redshift range of SNe la and extend this result to very high redshift GRB {z > 1.4) in a 
completely cosmological model-independent way. This approach has been also studied in [43J. 

Some analysis have been made with the GCGM and the GRBs as distant markers [44]. In the 
reference [35] the authors build a specific distribution of GRB to probe the flat GCGM and the 
XCDM model. While the GCGM has an equation of state given by expression (|2|) the XCDM 
model is considered in terms of a constant equation of state uj = p/p < 0. The main conclusion 
of this article is that the use of GRBs as a dark energy probe is more limited when compared to 
SNe la. We anticipate that we shall arrive at a similar conclusion. Moreover the XCDM model 
is better constrained than the GCGM. On the other hand, in [36] the GCGM and the ACDM 
model are compared by using the GRB and SNe la data to build the Hubble diagram. These 
authors show through the statistical analysis that the Chaplygin gas model (they use a = 1) 
have the best fit when compared with the data. Also they verify that the transition redshift 
between the decelerated and the accelerated state of the Universe occurs at z ~ 2.5 — 3.5 rather 
than z ~ 0.5 — 1 based on the analysis made with the SNe la. Here, for our purpose, we will 
assume the plausible assumption that GRBs are standard candles and we will use the data from 
Liang et al. [39] , calibrated cosmology-independently from the Union2 compilation of SNe la, to 
constraint the cosmological parameters of the GCGM. We want to show how GRB data could 
constraint different Chaplygin cosmologies. 

This paper is organized as follows. In next section, we described a brief review of GCGM. 
In section 3 the luminosity distance di is obtained for the GCGM and compared with the 
observational data. Finally, in section 4 we present our discussion and conclusions. 

2 The Generalized Chaplygin Gas Model 

We consider here an homogeneous and isotropic Universe described by the Friedmann's equation 



where the density p has the subscripts m for the matter pressureless fluid and c for the generalized 
Chaplygin gas with equations of state p m = and p c = —A/p®, respectively. Dot means 
derivative with respect to the cosmic time t. Flat, closed and open spatial sections correspond 
to k = 0, 1, —1 for the constant of the curvature. 




(3) 



3 



In our case, each fluid obeys separately the energy conservation law. The equations and the 
respective solutions are given by 



Pm + 3-pm = -> p m = , (4) 

a ( A\ ( - 1 - A \ 

A+3-(p c - — ) = -> / , c = Pc0 (^ + -^_ Y ) ( 5 ) 



where /3 m o = Pm(ao)> PcO = Pc(ao) = (A + B) l ^ 1+ °^ with a(i = 0) = ao = 1 being the scale 
factor today. The new definition of the constant A is given by A = Aj 'p^ a and it is connected 

to the sound velocity today in the gas by the expression v so = \J dp c /dp c = V aA. 

to 

Initially, the GCGM behaves like a dust fluid, with p oc a~ 3 , while at late times the GCGM 
behaves as a cosmological constant term, p oc Hence, the GCGM interpolates a matter 

dominated phase (where the formation of structure occurs) and a de Sitter phase. At the same 
time, the pressure is negative while the sound velocity is positive, avoiding instability problems 
at small scales [47]. 

In order to proceed with data comparison we need to calculate the luminosity distance in 
the GCGM. Using the expression for the propagation of light and the Friedmann's equation J3]), 
we can express the luminosity distance as 

d L = ^n = (l + z)S[f(z)} , (6) 
a 

where r\ is the co-moving coordinate of the source and 

S(x) = x for (Jfe = 0) , 
S(x) = sinx for (k = 1) , 

S(x) = sinhx for (k = —1) . (7) 

The function f{z) is given by 

1 f z dz' 

f{Z) = H~oJ {n m ( Z > + 1)3 + n c [A + (1 - A){Z' + l)3(l+a)]l/(l+a) _ ^ + 1)2} l/2 > ( 8 ) 



with the definitions 

_ 8^(2 p m0 _ 8ttG Pcq k_ 

m ~ Q tt2 i c ~ o tt2 ' k ~ tt2 ' v 9 / 



and Q m Q + O c o + fijt = 1. The final equations have been also expressed in terms of the redshift 
* = -l + 3- 

In our numerical calculations we relax the restriction that the pressureless matter component 
is entirely given by baryons. We consider the nucleosynthesis results for the baryonic component 
of the Universe and assume the total pressureless matter density as Q, m = Qb + Qdm, where 
n b h 2 = 0.0223 and H = lOOhKms- 1 Mpc' 1 . Then, in our notation f^m means the extra dark 
matter contribution. 



3 The Numerical Results 

The observational data set used in this article is composed by 42 GRBs from [39 \ 40]. As told in 
the introduction, this sample has been obtained and calibrated cosmology-independently from 
the Union2 compilation. This fact is of crucial importance to admit GRBs as cosmological 
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probes since the circularity problem described above is avoided. At the same time, this data 
set allow us to analyse the free parameters of the GCGM for a redshift range larger than the 
available data from SNIa reaching up to z ~ 6. It is important to emphasize that with this 
sample the authors of [39, 10 have obtained stark constraints on the Cardassian and Chaplygin 
gas model by combining the GRB data with other cosmological probes. 

If we want to have a reliable sample of GRBs to make our analysis, the Hubble diagram for 
the GRBs should be calibrated from the SN at z < 1.4. This allows to obtain the following 
luminosity/energy relations: The T[ ag — L relation, the V — L relation, the L — E p relation, the 
E~ — E p relation, the trt — L relation, the Ei so — E p relation, and the Ei so — E p — % relation. 
In general these relations can be written as log(y) = a + b log(x) (two-variable relations) and 
log(y) = a+bi log(xi)+62 log(x2) (multi-variable relation). In this relations y is the luminosity in 
units of erg s _1 or energy in units of erg and x is the GRB parameter measured in the rest frame; 
in the latter expression x\ and X2 are E p (l + z)/(300 keV) and + z)/(l day) respectively, 
and b\ and 62 are the slopes of x\ and X2 respectively. The calibration's process is achieved using 
two methods: linear interpolation (the bisector of the two ordinary least-squares) and the cubic 
interpolation (the multiple variable regression analysis). The variables a and bi are determinated 
with 1 — a uncertainties. With the linear interpolation, the error of the interpolated distance 
modulus can be calculated by a p = ([(z i+ i - z)/ (z i+1 - Zi)] 2 e 2 ^ + [(z - %)/ (z i+1 - Zi)] 2 e e ^ i+1 ) 1/2 , 
where \i is the interpolated distance modulus of a source at redshift z, i and e e ^ i+l are errors 
of the SNe, fn and are the distance moduli of the SNe at nearby redshifts Zi and Zj+i, 
respectively. In the case of the cubic interpolation method the error can be estimated by the 
expression <r M = (A^e 2 ^ + A 2 e 2 ^ i+1 + Ale 2 ^ i+2 + A^e 2 ^^) 1 / 2 , where e^i+j are errors of the SNe 
and are the distance moduli of the SNe at nearby redshifts Zi + j (index j run from to 3) 
with: 

[{Z i+ i - z)(z i+2 ~ z)(z i+3 - z)} 
[{Z i+ i - Zi)(z i+ 2 ~ Zi)(z i+3 - Zi)] 
[(Zj - z)(z i+ 2 ~ z)(Zj +3 - z)] 
[(Zi - Z i+ l)(z i+ 2 ~ Z i+1 )(z i+3 - Zi+l)] 

[{zj - z)(z i+1 - z)(z i+3 - z)] 

[(Zi - Zi +2 )(z i+ l - Zi +2 )(z i+3 - z i+2 )] 

[(Zj - Z)(z i+ l - z)(z i+ 2 ~ Z)] 
[{Zi - Z i+3 ){z i+1 - Z i+3 )(z i+ 2 ~ Z i+3 )} 

(10) 

The results obtained by the cubic interpolation method are almost similar to the results obtained 
by the linear interpolation method. It is important to emphasize again that the calibration 
results are completely independent of cosmological models used (for further discussion, see [38J ) . 

In order to compare the GCGM with the observational data, the first step is to compute the 
theoretical luminosity distance 

"'" = 51 ° g (^) +25 • <»> 

with the relations for the GCGM described above. Here, as in [391 140] ■ by using only linear 
interpolating we have the 27 GRBs at z < 1.4 from the Union SNe la data and the 42 GRBs at 
z > 1.4 obtained with the five relations (77^ — L,V — L, L — E p , E^ — E p , trt — L) calibrated 
with the sample at z < 1.4 that uses also the linear interpolation method. It is assumed 
that the GRBs luminosity relations do not evolve with redshift, so we could get the luminosity 
(L) or energy (-EL) of each burst at high redshift [z > 1.4). The weighted average distance 
modulus from the five relations for each GRB is \i = (J^ f^i / a fn) / a ^n)i with its uncertainty 
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Hi = c^ 2 ) -1 ^ 2 ) where the summations run from 1 to 5 over the five relations described 
above. 

Considering a set of free parameters {p} the agreement between theory and observation is 
measured by minimizing the quantity, 

^,^ MW! , (12) 

i=l a% 

where /j, th and fi obs are the theoretical value and the observed value of the luminosity distance 
for our model, respectively, and a means the error for each data point. We use Bayesian analysis 
to obtain the parameters estimations through the probability distribution function (PDF) 

P = Be~^, (13) 

where B is a normalization constant. A full Bayesian analysis is made by considering all free 
parameters of the model. However, we will study some particular Chaplygin configurations 
before a detailed analysis with 5 free parameters. With this strategy we hope to gain some 
intuition about the GRB data from the partial outcomes. Below, we will describe different 
Chaplygin-based cosmologies investigated in the present work. We show our results in Table 1-2 
and in Figures 1-7. 

Our first step is to study the Chaplygin gas (a = 1). We remenber that this equation of 
state, as cited above, has also raised interest in particle physics thanks to its connection with 
string theory and its supersymmetric extension |18j . We shall consider the prior information: 
< A < 1, < Q dm < 0.957 and < H < 100. The curves in Fig. Q] represent the 
99.73%, 95.45% and 68.27% contours of maximum likelihood after the first marginalization, 
i.e, integration of the likelihood function over the non requerid parameter. However, the most 
robust parameter estimation is the central value for the parameter obtained in the maximum 
of the one-dimensional PDF as in Fig. [2j From Fig. [2] we can obtain the final parameter 
estimation. We found Q dm = 0.04t^,A = 0.96±g;gf and H = 51.3±§;| at la level. However, 
the dispersion in the GRB data is quite high. We compute these same estimatives using the 
Supernovae Constitution sample [48J in order to compare the dispersion of these two samples. 
For the SN we found VL dm = 0.00±ftj£, A = 0.99±g;2i and H o = 59.7l 2 ^ at la level. Some 
constraints on the Generalized Chaplygin gas have been placed using the Constitution data set 
|49j . This allows a comparison between some of our results and the ones from Supernovae. In 
general, GRBs recover the results from SNe but with a high dispersion. 

In our next analysis we relax the prior information about the Hubble parameter and leave 
it free to vary. We show the two-dimensional PDFs in the Fig. ([5]). In Fig. ([6]) the solid lines 
are the corresponding one-dimensional probabilities. 

The above choice for the priors in the parameter a is conservative. With this choice we want 
to avoid a super luminal propagation in the sound speed formula. However, as argued in [50] 
the formula Vg = aA represents the group sound velocity. Actually, in order to violate causality 
the wavefront velocity should exceed 1 [51]. Considering this possibility we assume now a > 
and compute the one-dimensional PDF as showed in dashed lines in Fig. [6l 

Until now we have considered a flat Universe in our analysis. In order to have a more 
general statistical analysis we allow a non-vanishing curvature in our model. The complete five- 
dimensional analysis is computationally hard but still feasible. We assume as prior information 
that our Universe deviates slightly from the flat model assuming 0^ to vary between [-0.6,0.6]. 
For this case, we show the results in Fig. [71 
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Table 1: For the different Chaplygin-based cosmologies in the first column we show the final ID 
estimation for the free parameters. The errors are computed at la level. 
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Table 2: Best fit values for the 2D PDFs. 
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Figure 1: Two-dimensional probability distribution function (PDF) for the free parameters in 
the CGM. The curves represent 99.73%, 95.45% and 68.27% contours of maximum likelihood. 
The darker the region, the smaller the probability. 
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Figure 2: One-dimensional PDFs for the three free parameters of the CGM. 
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Figure 3: Two-dimensional PDFs for the GCGM fixing Hq = 72 km s Mpc . The curves 
represent 99.73%, 95.45% and 68.27% contours of maximum likelihood. The darker the region, 
the smaller the probability. 




Figure 4: One-dimensional PDFs for the three free parameters of the GCGM when Hq = 
72 km s~ 1 Mpc^ 1 . 



4 Discussion and Conclusions 

In this study we have analyzed the Chaplygin gas model with a sample of 42 GRBs. Although 
the use of GRBs as a cosmological tool is a promising way to probe cosmology at high redshifts 
we have verified that the available data is still insufficient to impose precise constraints in 
cosmological models. As observed in our analysis, the dispersion is still high when compared 
with others observational data sets. However we hope that with the future data from the 
final Swift BAT Catalog we will be able to put strong constraints on the dark energy/matter 
properties. 

In our analysis, the unification scenario was not imposed from the beginning. This means 
that we allow an extra dark matter contribution (f^m) i n our calculations in order to probe 
whether the unification scenario is favoured. In our first analysis the free parameters (A, Sldm 
and Hq) of the Chaplygin gas (a = 1) were well constrained. Our results are in agreement with 
the Supernova results [52]. The only difference is that we find a lower value for the Hubble 
parameter, Hq = 51.3jVy (lc). However, it is possible to find in the literature similar results 
for the parameter Hq [53J. 

In our second analysis, in order to check the behaviour of the model when Hq = 72 km s~ 1 Mpc 
we leave a free, that is the so called Generalized Chaplygin Gas Model. From Figs. [3] and U the 
unification scenario is again favoured. However, the uncertainties are still high. The parameter 
a assumes a large negative value. There is no any peak in the parameter a distribution and 



8 




40 50 60 70 35 40 45 50 55 60 65 70 35 40 45 50 55 60 65 70 



Hq Ho Hq 



Figure 5: The same as Fig. (1) but considering four free parameters for the GCGM and the 
prior < a < 1. 




Figure 6: One-dimensional PDF for the GCGM free parameters when Hq is free to vary. The 
solid lines correspond to the prior choice < a < 1 while dashed lines correspond to the prior 
a > 0. The final estimation for the parameter a does not depend on its prior information. 

the probability remains constant for negative values. For the background dynamics the region 
(a < —1) represents a behavior different from the matter dominated phase when structures start 
to form. On the other hand, negative values for a imply an imaginary sound velocity, leading 
to small scale instabilities at the perturbative level. Rigourously, the general situation is more 
complex: such instabilities for fluids with negative pressure may disappear if the hydrodynami- 
cal approach is replaced by a more fundamental description using, e.g., scalar fields. However, 
this is not true for the Chaplygin gas: even in a fundamental approach, using for example the 
Born-Infeld action, the sound speed may be negative if a < 0. Perhaps the restriction a > 
must be imposed for all observational tests. 

We work also with a set of four free parameters. Varying all four parameters, the preceding 
results are confirmed. Leaving the parameter Hq free to vary, we confirm that the hypersurface 
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Figure 7: One-dimensional PDF for the GCGM free parameters when we allow a non-vanishing 
curvature. The solid lines correspond to the prior choice < a < 1 while dashed lines correspond 
to the prior a > 0. The final estimation for the parameter a does not depend on its prior 
information. 

Hq = 72 km s _1 Mpc -1 doesn't represent the maximum probability in the 4-D parameters phase 
space. Chaplygin gas models show values lower than Hq = 72 km s~ 1 Mpc~ 1 [52J. We observe 
also that there is a significant difference in the final parameter estimation when we consider the 
prior a > 0, instead of < a < 1. For instance, the unification scenario (fldm = 0) is favoured 
only with the choice a > 0. Moreover, there is now a peak in the a distribution at a = — A.3^f 2 
but with a high dispersion. Again, negative values for a are favored, despite the two-dimensional 
PDF (a x Hq) in Fig J5] indicates a high probability for a > 6. Such a contradiction seems to be 
an artifact of the marginalization process as can be seen also in the two-dimensional PDFs (£ldm 
x Hq) and (A x i^o) m Fig. [5j These plots confirm that the final ID estimation can be very 
different from the partial 2D ones. This diferrence is due to the integration of the probability 
function over the adopted prior values of the remmaning parameters. 

The analysis with five free parameters confirm some of the previous results. Negative cur- 
vature is prefered as well as Sn data data [52]. Also, the parameter a is now estimated with a 
positive value, in constrast with the previous results. 

The Chaplygin gas parameters have been estimated in many papers, considering different 
analysis and several observational data sets. Constraints critically depend on whether one treats 
the Chaplygin gas as true quartessence (replacing both dark matter and dark energy) or if one 
allows it to coexist with a normal dark matter component. The former situation is widely 
considered in the literature. As we leave the density parameter Qdm free to vary in all the 
cases analysed here, it is not possible to directly compare our results with unified Chaplygin 
cosmologies unless we assume the prior Qdm = 0. This case has been studied using GRBs and 
other probes in reference [ID]. For a comparison with this reference, figure [8] shows the two- 
dimmensional probability for the free parameter of the unified (CI dm = 0) GCG model. The best 
fit occurs at (a = 0.15, A = 0.75). This result agrees (at la) with the joint analysis showed in 

The analysis of section 3 can be compared with [52], where the influence of a free Cldm 
parameter on the final estimations was taken into account. Our results have high confidence 
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Figure 8: Constrainst on the free parameters of the unified (Qdm = 0) GCG model. 



with the results obtained in [52J. 

Finally, we remark that, perturbative analysis of Chaplygin models, for instance, reveals 
a large positive value (a >> 200) for the parameter a |54] while kinematic tests show values 
negatives or close to zero. At the background level, the crossing of different data sets (including 
for example SNe, CMB, BAO, H(z) data and galaxy cluster mass fraction) will provide a more 
accurate scenario for each Chaplygin-based cosmology studied in this work. We leave this 
analysis, including the perturbative study, for a future work. 



H Marginalized 
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